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Abstract. The luminosity distance - redshift relation for a wide class of 
generalized Randall-Sundrum type II brane-world models with Weyl fluid is 
compared to the presently available supernova data. We find that there is a 
class of spacially flat models with different amounts of matter Qp and Weyl fluid 
fid, which have a very similar fitting quality. The best-fit models are equally likely 
and can be regarded as extensions of the ACDM model, which is also included. 
We examine three models with different evolutionary history of the Weyl fluid, 
characterized by a parameter a = 0, 2 and 3. The first model describes a brane 
which had radiated energy into the bulk some time ago, but in recent times this 
energy exchange has ceased and only a dark radiation (o = 0) is left. In the other 
two models the Weyl-fluid describes a radiating brane throughout the cosmological 
evolution, up to our days. We find that the trought of the fitting surface extends 
over a wider Q^-range with increasing a, but the linear correlation of and Qp 
holds all over the examined Qd range. 



1. Introduction 

Current observational data [l]-[3] suggest that the cosmological model of a Universe 
with only baryonic matter has to be modified. In the easiest way, the model can 
be reconciled with observations by the introduction of a cosmological constant A 
and of considerable amount of dark matter (ACDM model). Because the energy 
densities of both baryonic and dark matter decrease during cosmological evolution, 
the cosmological constant will dominate the late-time evolution. This process was first 
suggested for the explanation of la supernovae data, which suggest that our Universe 
has reached an accelerating phase. In a A-dominated universe, the luminosity distance 
increases faster with redshift than in the model without A [4j , exactly as required by 
the supernova data. 

Generally, the agreement with experiments can be achieved by introducing a dark 
energy component of the Universe, which replaces A. Such a dark energy in general 
does not clump. A recent analysis [Fj shows that a dark energy model with varying 
dark energy density going through a transition from an accelerating to a decelerating 
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phase at redshift 0.45 fits well the observational data. Based on observations, the dark 
energy equation of state w ^ p/p is within about — 1 ± 0.1 [16]. 

It has been expected for some time that alternative gravitational theories, 
motivated by string / M-theory could replace dark matter and dark energy by 
geometric effects. The curved generalizations (see for example the review [7]) of the 
original Randall-Sundrum type II model [8] consist of a hypersurface with tension A 
(the brane), representing our observable universe, embedded in a 5-dimensional space- 
time (the bulk). Gravitational dynamics on the brane is governed by an effective 
Einstein equation [S], [TU]. The sources of gravity in the effective Einstein equation 
include terms due to the asymmetric embedding of the brane into the bulk jlO l , non- 
standard model fields in the bulk, and even quantum corrections approximated as 
induced gravity effects [l2]-p!5]. 

The most relevant source term for early cosmology is a quadratic source term in 
the energy- momentum tensor [TT]. This term dominates over the linear term before 
the Big Bang Nucleosynthesis (BBN) . In the simplest case of cosmological symmetries 
and suppression of the energy exchange between the brane and the bulk and whenever 
the bulk contains a static black hole, the Weyl curvature of the bulk generates a 
so-called Weyl fluid effect on the brane. In Fig 1 of our companion paper [16] (to 
be referred in what follows as paper I) we classify the different brane-world theories 
and their inter-relations. They are divided into two branches, one containing the 
original Randall-Sundrum type II model (BRANEl) and the other the flat DGP model 
(BRANE2). The model with Weyl fluid belongs to the BRANE 2 branch. 

Supernova data were confronted with the induced gravity models [l3-[2D]- When 
they are combined with the Sloan Digital Sky Survey (SDSS) baryonic peak, these 
seem to rule out the flat DGP models [T^ , [H] • However it was argued in [^Oj that the 
Cosmic Microwave Background (CMB) shift parameter can over-turn this conclusion. 
Structure formation and CMB were also considered in the DGP models in Ref. [21]. 

Most recently, the authors of [22] tested the accelerating phase of the universe's 
expansion with a comparison of the models and the supernova data. They have tested 
the ACDM model, the DGP model and three wCDM models with equations of state 
where w{a) (i) was constant with scale factor a, (ii) varied as w{a) — wo + Wa{l — a) for 
redshifts probed by the supernovae but fixed at —1 for earlier epochs, and (iii) varied 
as Wo + Wa(l — o) since the recombination. Their main conclusion is that all the five 
examined models explain equally well the acceleration, and none of them could be 
selected as a preferred model, based on the la-type supernova data. 

The authors of Ref. [TB] have compared the predictions of the flat BRANEl 
and BRANE2 models to the Gold [23] and Supernova Legacy Survey (SNLS) [24] 
supernova data sets, incorporating the baryon acoustic peaks into the analysis. These 
brane-world models in certain parameter range (when their induced gravity parameter 
fli is small; the flat DGP models falling outside this range) are satisfied by both 
supernova data sets. The BRANEl models fit better to the SNLS data, while the 
BRANE2 models fit better to the Gold data set. Since the analysis depends very 
weakly on the bulk cosmological constant A, the value of A was fixed at zero. With 
this modification, the BRANEl model fits better to the SNLS data than the ACDM 
model and fits comparably well to the Gold data. The same conclusion holds for the 
BRANE2 model. In the analysis of [TH] the dark radiation dimensionless parameter 
Qd is switched off. 

Using two recent supernova data sets, the CMB shift parameter, and the baryon 
oscillation peaks, the authors of Ref. [25] have found that the LDGP model (a 
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subclass of the BRANEl models with the effective energy density having a phantom- 
like behavior due to extra-dimensional effects, see Fig 1 of Paper 1) fits the observations 
if it is very close to the ACDM model. The modification of the LDGP model with 
respect to the ACDM model appears in the form of a linear term in the Friedmann 
equation, H/rc, where H is the Hubble parameter and Tc a crossover scale. This 
model includes a cosmological constant, possibly screened by the modified gravity, 
however the comparison with observations sets strong constraints on the screening. 

The first comprehensive study of the generalized Randall-Sundrum type II (RS) 
brane- worlds tested against astronomical data was presented in Ref . [55] . Agreement 
with earlier supernova data has been established in the presence of a cosmological 
constant. In this analysis the dark radiation from the bulk was switched off (fid = 0) 
and the energy-momentum squared term was kept. Under these assumptions, for flat 
spatial sections and matter parameter Qp — 0.3 the maximum likelihood method gave 
fl\ = 0.004 ± 0.016 for the parameter characterizing the source term quadratic in the 
energy-momentum. This in turn implies a tiny value of the brane tension, which is 
disfavored by generic brane-world arguments. Moreover, much lower values for ^\ 
emerge from both CMB and BBN. 

In contrast to Refs. [TS] and the analysis of [57] keeps both and ^l^, 
the latter obeying \Qd\ < 0.01. The best fit is obtained at Hp = 0.15, Ha — 0.80, 
r^A ~ 0.026 and Hd — 0.008. With the high value of the brane tension set by either 
(a) the value of the 4-dimensional Planck constant and sub-millimeter tests [21] on 
possible deviations from Newton's law (in units c = 1 = h these give ^'^^^letop — 138.59 
TeV' see [29], [7]), (b) astrophysical considerations X'^Zo = 5 x 10® MeV^ ^ or (c) 
BBN constraints A^^^ — 1 MeV''' [31], the quadratic source term barely counts at 
late-times in the cosmological evolution. 

Given the high limits for the values of A, in any realistic model Qx can be safely 
ignored. This is a crucial difference of our forthcoming analysis as compared to the 
one presented in Refs. [27] and [26], where the corresponding cosmological parameter 
fix was kept. 

The next question is whether the source term arising from the Weyl curvature of 
the bulk may be kept, in other worlds, whether Qd 7^ 0. The Weyl curvature of the bulk 
gives an energy density pd = 6m/K^a'*, where = SttG is the gravitational coupling 
constant. In Ref. [32] it was shown that the BBN limits constrained the dark radiation 
component as —1.23 < Pd (zbbn) / Pj (zbbn) < 0.11. Combining this with CMB 
constraints reduces this range to —0.41 < pd {zbbn) I Pi (zbbn) < 0.105. Here p^ is 
the energy density of the background photons. Another constraint for the value of the 
dark radiation at BBN was derived in [S^ as —1 < pd {zbbn) / Pv {zbbn) < 0.5, where 
pi, is the energy density contributed by a single, two-component massless neutrino. 
This constraint was derived for high values of the 5-dimensional Plank mass. 

In the simplest case the Weyl source term evolves as a radiation, thus its present 
value is obviously tiny. This is the reason why all mentioned references [18] and 
|26j comparing RS brane-worlds with observations disregard dark radiation. But is 
this a necessary assumption? Formulating the question the other way around: if we 
include even a small component of dark radiation into the late-time universe model 
we face a serious problem. Due to the fact that the energy density of dark radiation 
decreases as a^"* (compared to that of matter which is a~'^), even an amount of 
dark radiation of the same order as the amount of baryonic matter nowadays implies 
dark radiation dominance in the past, for example during structure formation. This 
conclusion is contradicted by numerical simulations, which favorize cold dark matter 
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as the dominant component of the Universe during structure formation [34] . 

However we can generaUze the vaUdity of the model by Hfting the requirement of 
a constant mass m in the dark radiation energy density. A constant m imphes a static 
Schwarzschild-anti de Sitter bulk and no energy exchange between the brane and the 
bulk. Therefore dark radiation is a manifestation of an equilibrium configuration with 
a static bulk, and it may be well possible that such a situation is reached only at the 
latest stages of the evolution of the brane-world Universe. Whenever m depends on a 
certain, non-zero power of a, the evolution of the energy density of the Weyl source 
term evolves in a non-standard way, allowing to escape from the argument of a small 
Weyl fluid left nowadays. 

We propose here the LWRS (Lambda-Weyl fluid-Randall-Sundrum) model, a 
specific RS model with i) cosmological constant, ii) the brane radiating away energy 
during various stages of the cosmological evolution, characterized by the index a 
and Hi) a Weyl fluid depending on the actiul value of a during the latest stage 
of cosmological evolution^ which can be tested by supernova observations. For the 
inclusion of the LWRS model in the classification of brane-world models, see Fig 1 of 
paper I. 

The LWRS model takes into account the possibility of an energy exchange between 
the brane and the bulk. This idea is not new. Indeed, it was already proposed 
that during an inflationary phase on the brane radiation is emitted and black holes 
thermally nucleate in the bulk [35] . Later on, but still in the high energy regime, the 
brane radiates such that the mass function of the bulk black hole increases with 
[36] . This means that the Weyl source term becomes a constant in this era. The brane 
continues to radiate away energy during structure formation [37] , a process leading to 
a bulk black hole mass function m (x a", with 1 < a < 4. (Other models with the 
brane radiating energy into the bulk are also known [38".) 

For a = the Weyl fluid is known as dark radiation, for a = 2, 3 it gives the 
correct growth factor during structure formation. For a — 1, 4 it is indistinguishable 
from dark matter and a cosmological constant, respectively. Therefore pure dark 
radiation can emerge only in the low-z limit, while at earlier times a dynamic bulk - 
brane interaction governed by energy exchange should be present. 

2. Confronting the models with the selected supernova data 

As type la supernovae result from the explosion of white dwarf stars with identical 
mass, they show remarkable similarities. By employing well established calibration 
methods, one can calculate the maximal luminosity of the object (in the reference 
system of the explosion) . This is done by analyzing the time-dependent variation of 
the emitted luminosity and the spectrum, a method known as the Multi-Color Light 
Curve analysis [39] . [23| . In this process the observed parameters, the shape of the 
light curve and the spectral distribution of the emission have to be converted into 
the reference system of the host galaxy. For distant supernovae this translates to 
take into account the time dilation and the so-called i^-correction [40]. While these 
methods depend on z, they are independent on the specific cosmological model. After 
performing these corrections, we have well-calibrated maximal luminosities for the 
supernovae of type la and in consequence they are considered as standard candles. 

In 2003 a hst of dL-z data pairs were pubhshed for 230 supernovae of type la [IT] , 
and 60 of them had low absoprtion (i.e. Ay-l) and cosmological redshift (z > 0.01). To 
give result which can easily be compared to earlier works, we also involve this selected 
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low-absorption data set for the most of the examinations. The basic supernova data 
we use here is the improved Gold set |35], which was released in 2006. 

We confront with supernova observations several models from paper I. In Fig[T]we 
represent graphically on both logarithmic and linear scales their luminosity distance - 
redshift relations up to z = 2.5 . The plots are for fc = and flp = 0.27 (according to 
the combined analysis of the SDSS and WMAP 1-year data in Ref 2 ). In particular, 
the luminosity distance - redshift relation is shown for the following models: 

• The LWRS model (the perturbativc solution given by Eqs. (56), (60)-(61), (63), 
(65) and (67) of paper I, with Qx = and for a = 0) for the two values of 
the late-time dark radiation Qj. = —0.05 and = 0.05 (the curves 1 and 3, 
respectively). The latter models contain a brane which radiates energy at early 
times (for f2d > 0) and during structure formation, such that a bulk black hole 
is formed and its mass increases continuously. As this process slows down, the 
Weyl curvature of the bulk induces the late-time dark radiation on the brane. 

• The ACDM model, given by Eqs. (60)-(61) of paper I (curve 2). 

• The solutions with brane tension A = 2A/k^ and no dark radiation (given by Eq. 
(52) of paper I) for both admissible values for this model, at JIa — 0.704 (curve 4) 
and f^A = 0.026 (curve 6). The former is similar to the class of models discussed 
in [26] . 

• The late-time universe fl\ = limit of the RS model with Randall-Sundrum fine- 
tuning, containing a huge amount of dark radiation ild = 0.73, given by Eq. (44) 
of paper I (curve 5). 

In Fig. 1 we plot these models in a comparison to low-absorption supernova 
data from Ref. [41] (red triangles) together with the Gold set [23] (black dots). The 
error bars are indicated in the respective colors. The diagrams with linear scale are 
more instructive, as they emphasize the difference among the predictions of the chosen 
models and how they fit data, while the logarithmic scale better disseminate between 
the low z points. 

The models represented by the curves 1, 3 and 4 by eye seem to compare as well 
with the supernova observations as the ACDM model (curve 2). By contrast, the 
models represented by the curves 5 and 6 seem to be not supported by observations. 
The model with no cosmological constant and significant dark radiation fl^ = 0.73, 
r^A = (curve 6) and the model with A = and = 0.025 (curve 5) are 

significantly inconsistent with the observations, as they give = 213 and 395, 
respectiveljQ. All other models shown on Fig [T] are comparable with the supernova 
observations, as it was expected by a simple glance. 

The — 50 value found for the A — K^A/2-model with f^A — 0.74 (curve 4) is 
slightly better than found the ACDM model. However, as mentioned earlier, the 
tiny brane tension A — 38.375 x lO^^^TeV*, several order of magnitudes lower than 
all existing lower limits rules out this model as well. 

The best fitting models are the models with brane cosmological constant; a high 
value of the brane tension (leading to fl\ « 0) and a small contribution of dark 
radiation, fid = ±0.05 (the curves 1 and 3). For ild — —0.05 we find ~ 65, which 
is still acceptable. For fid — 0.05 we get = 49. 

I We mention here that we have also excluded several other models with Randall-Sundrum fine- 
tuning (not shown on Fig[T]l, which have either a very low value of the brane tension or a significant 
dark radiation. For example, the models with = 0.0258835, Qx = 0.70412 and = 0.70412, 
Ux = 0.0258835 gave = 246 and 415, respectively. 
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Figure 1. (Color online) Luminosity distance - redshift relations for selected 
brane-world cosmologies and for the ACDM model, compared to the supernova 
data. The diagrams are log-scaled (left panel) and linearly scaled (right panel). 
Selected low absorption supernovae from Ref. 1411 are plotted with red, black dots 
represent the Gold set '>'■'> . Both sets are represented with the corresponding error 
bars on the log-scaled diagrams. For the sake of perspicuity, the error bars of low 
absorption supernovae are not represented on the linearly scaled diagrams. The 
plotted models are the ACDM model (2); the brane models with cosmological 
constant and late-time dark radiation (1 and 3); without cosmological constant 
but with dark radiation (5); with cosmological constant satisfying A = 
thus low brane tension (4 and 6) and no dark radiation. 



Values of fid between these limits are also admissible. It is likely that by increasing 
fid towards higher positive values, remains compatible, however the accuracy of 
the perturbative solution is deteriorated with increasing fid, therefore higher orders 
in the expansion would be necessary to take into account. 

3. The Gold2006 set of supernovae 

More recently, Riess et al. [42] have published a new set of 182 gold supernovae, 
including new HST observations and recalibrations of the previous measurements. It 
is an interesting question how this rccalibration influenced the above conclusions for 
the well-fitting models with dark radiation. 

We applied the same tests to the Gold2006 data set as described in the previous 
section. First we assumed that flp — 0.27, as before, cf. Ref [2]. In this case the 
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Figure 2. (Color online) The luminosity distance - redshift relation for the viable 
branc- world models and the ACDM model (the curves (l)-(4) of Fig[TJ, both with 
logarithmic (left panel) and linear scale (right panel), compared to the smeared 
Gold set 1231 . The best fit is obtained for the brane- world model (3), with 5% 
dark radiation. 



critical value of is 197 at 80% level and 209 at 90% confidence level. Then the 
models represented by the curves 1-4 of Fig [2] behave as follows. The model with a 
small amount of negative dark radiation is disfavored at 80% confidence (x^ — 204). 
The models with A = 2K/ and J^a — 0.704 are ruled out at 90% confidence level, too 
(as x^ = 221). As expected from the previous analysis, the ACDM model (x^ = 192) 
and the LWRS model with il^ = 0.05 (giving x^ — 194) compete closely. We also 
remark that varying 0,^^ between —0.03 and 0.07, the x^ remains under the critical 
value. 

For gaining a deeper insight we have then calculated the predictions of the models 
between Vld = —0.10 0.10 with a stepsize of 0.01 in fi^, with Vlp allowed to freely 
vary in the domain 0.15 0.35 and z in the range 0-^3. Then we looked for the best 
fit of the Gold2006 set in the Vld — Vlp space. This is represented on FigO The global 
minimum of the surface is at Q.d — 0.040, ilp = 0.225 (x^ — 190.52), which suggests 
an interesting opportunity for a Universe with less baryonic density and with dark 
radiation, compatible with the Gold2006 supernova data. The 1-a confidence interval 
is centered about this value. The ACDM model (where Q.d is exactly 0) has the local 
minimum of Vlp = 0.275 (x^ — 195.8), but this is outside the 1-a confidence interval. 

Similar conclusions emerge from the plot in the n\ — flp plane. Fig U) Here the 
global minimum of the surface is at JIa = 0.735, flp = 0.225. The local minimum of 
the ACDM model is at = 0.725. 

We note that there is a forbidden parameter range in both planes fid — and 
Q\ — Qp, represented by white regions on Figs [3] and H) This is because the Friedmann 
equation for these brane-world models 



H{z) 



combined with fl\ 



Ho 



f^A + np{i + zf + nd{i + zf > 



Qd = 1 gives the constraints 



"(l + ^)'-l" 
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1 > 



(1) 
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Figure 3. (Color online) The fit of the luminosity distance - redshift relation for 
the LWRS brane-world models with dark radiation, (including the ACDM model 
for = 0). There is no assumption for Qp £ (0.15, 0.35), its preferred value 
0.225 being determined from the supernova data, together with the preferred 
value 0.040 of Qj. The contours refer to the l-a and 2-a confidence levels and 
both are centered on the LWRS model with the values given above. The local 
minimum represented by the ACDM model is at flp = 0.275. Both the global and 
local minima are marked. The white area in the lower left corner represents the 
forbidden region of the parameter space for 2 = 3. 



ii + zy[i + zii-np)] <o 



in the Qd — plane and 

in the Qa — plane. 

The forbidden region increases in both cases with z. 
the limits to z cx), we obtain the limiting curves limz 



(3) 



If we would like to extend 
^oo f^rf'" (z, ^p) = in the 
Op plane and limz^oo f^™'*'^ {z, 17p) = 1 — fip in the J^a — ^p plane. However the 
LWRS model being valid only for low values of z, we represent on the graphs only the 
forbidden range for z = 3. 



4. The compatibility of the LWRS model 17^ 
cosmological evolution 



0.04 and a = with 



The energy density of dark radiation decreases too fast during cosmological evolution 
to result in a considerable amount nowadays. We discuss this problem and its possible 
remedy in detail here. First we comment on the problem, then we show how an energy 
exchange between the brane and the bulk can leave a considerable amount of dark 
radiation. 
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Figure 4. (Color online) Same as on Fig|3l but in the Q\ — Qp plane. The global 
minimum is at Q\ = 0.735, Qp = 0.225, while the local minimum for the ACDM 
model gives f2yv = 0.725 and Qrn = 0.275 (both marked). The white area on the 
top right corner represents the forbidden parameter range. 



The constraint derived in [32] for the energy density of the dark radiation: 

_ 0.41 < ^4^^ < 0.105, (4) 

[zbbn) 

where p-y (zbbn) — PT^^j^ is the energy density of the background photons at the 
beginning of BBN. The coefficient 



'^B 



3.78 X 10"^^ J m"-^ 



(5) 



30^^(fc)3 

contains [11], [43] the effective number t/, of relativistic degrees of freedom, which 
depends on the temperature. According to [44^ = 10.75 at the beginning of BBN, 
when Tbbn = 1-16 x 10^° K. Thus p-y {zbbn) = 7.37 x lO^'' J m'^ emerges, giving 
the constraint 

- 3.02 X 10^5 Jm~3 ^ (zbbn) < 7.74 x 10^^ Jni"3_ 

Note, that the domain of aUowable negative values is larger than the one for positive 
values. 

As for today the background photons have cooled to Tq = 2.725 K and for such 
low temperatures = 3.36 [43], [44] their energy density p^ — p^ {z — 0) is 

p^ = 7.01 X 10"" Jm"^ . (7) 

With the value Hq = 73^1 km s^^ Mpc^^ of the Hubble constant J\, cf. Eq. (23) of 
paper I the present day cosmological parameters p and Q (both for background and 
dark radiation) relate as 

Pd,7 = 9.00 X 10~^° r2d,7 Jm"^ . (8) 
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Thus the present value of il^ is 

Q-y = 7.74 X 10"^ , (9) 

which is quite neghgible. If the Weyl source term were to evolve as radiation, its value 
would be even smaller, cf. Eq. (jl]). Indeed Eqs. ^ and ^ imply 

- 1.02 X 10"'' <fld< 2.62 X 10"^ . (10) 

\fld\ is of the same order of magnitude or smaller as Q^. 

However if the brane is radiating during structure formation, the mass parameter 
m becomes a function of the scale factor m oc a", with 1 < a < 4 [37]. Then the 
energy density scales as a**"". 

Now let us suppose that the brane is in an equilibrium (non-radiating) 
configuration with a = in the domain < z < 2:1. In a preceding era zi < z < 
the brane radiates such that a ^ 0, finally right after the beginning of BBN, 
at < z < zbbn there is equilibrium once more {a = 0). Here zbbn = 
{Tbbn/Tq) - 1 = 4.26 X 10^. According to this evolution 

' \ 4 / \ 4— Q / 

, , 1 oo \ / oi \ / a* 

Pd (Zbbn) — Pd \ — 



ai J \a*/ \aBB 



N 

Inserting this in Eq. ^ and employing Eq. ^ we obtain: 

- 1.02 X 10"'^ < f j^^^ ] < 2. 62 X 10"^ . (12) 



In the particular case a = we recover the constraint pop set on pure dark radiation. 
However for any a > we get 

2* > (1 + zi) [max (-0.98 Qd, 3. 82 1}^)]^^" x 10*^" - 1 . (13) 

Let us specify this result for the best fit value fid — 0.04. Depending on a 
we obtain the following numerical relations between the redshifts characterizing the 
switching on and off of the radiation leaving the brane: 

1527.80+ 1528.80 21 , a = l 

. 38.10 + 39. 10 zi , a = 2 

- ^ 10.52+ 11.52 zi , a = 3 ■ ^ ' 

5.25 + 6.25 ^1 , a = 4 

It is evident that the value of increases with zi (this dependence becoming an 
approximate scaling for higher values of zi) and decreases with a. 
The lower limit in the LWRS model is zi =3. Then 

6114.20 , a = l 

^. 155.40 , a = 2 

45.08 , a = 3 ■ ^^^^ 

24.01 , a = 4 

For the higher values of a the duration of the radiative brane regime necessary to 
produce a high value of fid today is quite short. 
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Figure 5. (Color online) Same as on Fig[3] but for the a = 2 models. 

5. LWRS models with a = 2, 3 confronted with supernova data 

In the absence of a known mechanism for changing a, we examine here the cases when 
a = 2 and a = 3 hold throughout the cosmological evolution, up to nowadays. For 
this we confront these models with the Gold 2006 set exactly as described before. To 
preserve the validity of the perturbative solution, the range of fid was selected to be 
—0.1-0.1, and we probed the range 0.15-0.35 of ftp. The assumption for flatness was 
kept, too. 

The results are qualitatively similar to the a = case. The remarkable difference 
is that the peak of the minimum turned into a "trough" , which lies aslope in the flp-fld 
space. This means that instead of a district solution, a complete model family exists 
in both cases, which can equally well explain the supernova data. The flp dependence 
of fid is less in the a = 2 model as compared to a = 0, and is very small if a = 3. The 
fid — case is the ACDM model where these model intersect. The steeper slope of the 
minimum trough thus allows a much lower range of flp in the a = 2, and especially 
in the a = 3 models, with a value close to 0.3. On the other hand, the range of fid 
gets more and more wide with increasing a, which results in the conclusion that the 
presence of fid is mathematically plausible, and they have to be accounted for in RS 
cosmology. 

Due to the higher slopes of the l-u and 2-a contours in these a — 2, 3 models, 
flp is much less affected by the Weyl fluid, while fid can have various values in the 
detriment of fl\. Therefore the Weyl fluid can explain some of the dark energy. 

6. Conclusions 

The luminosity distance given in paper I as function of redshift in terms of elementary 
functions and elliptical integrals of first and second type for various brane- world models 
with Weyl fluid was confronted with the available supernova data sets, including the 
Gold2006 data [42]. The tested models were: 
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Figure 6. (Color online) Same as on Fig[3] but for the a = 3 models. 

(A) The models with Randall-Sundrum fine-tuning, discussed in section 4 of paper 
I, with a considerable amount of dark radiation as a bulk effect, and a high value of 
the brane tension. 

(B) The two models discussed in subsection 5.1 of paper 1, which obey A = 
have no dark radiation and were integrable in terms of elementary functions. 

(C) The LWRS models (subsection 5.2 of paper I), with a brane cosmological 
constant, for which the luminosity distance could be given analytically as function of 
redshift to first order accuracy in the dark radiation. (Due to its smallness, the source 
term fix quadratic in the energy density was suppressed in the perturbative models 
of paper I.) 

The brane-world models (A) although interesting for historical reasons, do not 
comply with observations. Even if we introduce an extremely high amount of dark 
radiation 17^ = 0.73, tentatively replacing the cosmological constant in the energy 
balance Q\ + Qp + Qd + ^\ = 1, these models are quickly outruled by supernova data 
(curve 5 of Fig[T]). Dark radiation is not capable to replace the cosmological constant 
in producing a late-time acceleration, since it scales as usual radiation. The more 
we go back in the past, the higher becomes its domination over matter. Therefore 
a cosmological constant or dark energy is still needed in the generalized Randall- 
Sundrum type II models. 

Our analysis has also dismissed immediately the model (B) with fljy = 0.026. 
Surprisingly, the other toy model (B) with J^a — 0.704 was in good agreement with 
the Gold2006 data, but ruled out by its low value of the brane tension, similarly as 
the models discussed in Ref. j26j . A low brane tension is in disagreement with various 
upper limits set by cosmological and astrophysical tests. 

The perturbative approach of subsection 5.2 of paper I can be considered valid 
for a Weyl fluid with —0.1 < fid < 0.1. In this range the LWRS brane-world models 
(C) were confronted with supernova data and for a = the dark radiation with 
significant negative energy density ruled out. The fact that a positive dark radiation 
(corresponding to a bulk black hole rather than to a bulk naked singularity) is favoured 
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by the presently available best supernova data is in accordance with the early behavior 
of the RS model with late-time dark radiation, where the brane radiating away energy 
in early times leads to a black hole, which can further grow during structure formation. 

The remaining LWRS brane- world models with a = and Qd between —0.03 and 
0.07 (and changed accordingly) turned out to be excellent candidates for describing 
our universe, as they show remarkable agreement with the Gold2006 supernova data 
sets. If f2p is allowed to vary in the range (0.15, 0.35), the preferred values are 
fid = 0.040, fip = 0.225, flA = 0.735. 

The preferred cosmological parameters determined by comparing the LWRS 
model with a = with supernova data alone are in perfect accordance with the WMAP 
3-year data. Indeed according to Ref. [3] ilph^ = 0.127toml and h = 0.73^55;^^ from 
which Qp = 0.238l^;^!^^ emerge. The value of flp determined by comparing the LWRS 
model with the supernova data alone is well in the middle of the domain allowed by 
the WMAP 3 year data. 

We have then proved that the preferred value of fid = 0.04 is compatible with the 
known history of the Universe if the brane radiates away energy into the bulk during a 
relatively short period of the cosmological evolution. Such a process occurring between 
2 = 24 and z — 3 could increase the amount of dark energy today with a factor of 
10'^ as compared to the non-radiating brane, exactly as required by the LWRS model 
with a ^ 0. 

The LWRS models with a = 1 (a = 4) are identical with the ACDM model 
with the only difference that some fraction of the dark matter (of the cosmological 
constant) has geometric origin. 

Finally, the LWRS models with a = 2 and a = 3 do not present a sharp 
minimum, but rather an elongated trought shape in the parameter space, with the 
slope increasing with the value of a and flp « 0.3. This means that in this class of 
models a wide range of values for fid (with a slight preference for negative values) and 
corresponding values for Q\ are fitting to the supernova data. 

We must note that the reliability of these values is somehow deteriorated by the 
relatively small number of high-z supernova and by the inherent difficulties in the 
calibration of the available data. An obvious source of error is that data from the 
Gold2006 set is a combination of measurements taken on different instruments [IS] 
and in fact it has been already signaled that the Gold2006 data set is not statistically 
homogeneous j46] . 

The conclusion of this paper is somewhat similar to that of Ref. [22] : the presently 
available supernova data are not enough to discern among several cosmological models. 
However the difference between the predictions of the acceptable models of our analysis 
(the ACDM model, the LWRS brane-world with a = and fid = 0.04 and the models 
with Weyl fluid and a = 2, 3) are increasing with z. One may reasonably hope that 
the very far (z > 2) supernovae, which will be discovered for sure in the following 
decade, will improve their comparison. 
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